Space-adiabatic Decoupling to All Orders 
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A systematic perturbation scheme is developed for approximate solutions to the time-dependent 
Schrodinger equation with a space-adiabatic Hamiltonian. For a particular isolated energy band, 
the basic approach is to separate kinematics from dynamics. The kinematics is defined through a 
subspace of the full Hilbert space for which transitions to other band subspaces are suppressed to 
all orders and the dynamics operates in that subspace in terms of an effective intraband Hamilto- 
nian. As novel applications we discuss the Born-Oppenheimer theory to second order and derive 
the nonperturbative definition of the g-factor of the electron within nonrelativistic quantum elec- 
trodynamics. 

PACS numbers: 03.65.-w, 03.70.-(-k 



The importance of slow parameter variations is insep- 
arably linked to early quantum mechanics |Q] and stood 
godfather for the time-adiabatic theorem of Born and 
Fock 1^ , later rigorized by Kato . The physical idea is 
simple and best demonstrated in the context of the dy- 
namics of molecules. The nuclei move slowly and follow, 
within a good approximation, some classical trajectory. 
The electrons are then subject to a potential of slow time 
variation. Abstracting from the particular setting, one 
can think of a time-dependent Hamiltonian H{t) acting 
on the Hilbert space TL. The slow time-dependence is in- 
troduced through H{st), where £ <SC 1 is a dimensionless 
scale parameter, and one wants to study the solution to 
the time-dependent Schrodinger equation 



at 



(1) 



in the limit of small s. For simplicity ft = 1, except 
for the physical examples. In (|^) we followed the stan- 
dard practice to switch to the slow time scale through 
the substitution t/e for t. It should only be recalled that 
t of order 1 corresponds to very long times when mea- 
sured in atomic units. We assume that H{t) has a phys- 
ically distinguished "relevant" subspace of dimension m 
with, only for the sake of presentation, single energy E{t), 
which is isolated from the rest of the spectrum. The pro- 
jection P(t) onto the relevant subspace is spanned by 
the instantaneous eigenvectors \(pa{t)) oi H{t) satisfying 
H{t)\(pa{t)) = E{t)\ipa{t)), a = 1, m. The initial wave 
function, ipo, is assumed to lie in the relevant subspace, 
■00 = PiQ)4'o- There is no reason to expect persistence 
as ipt = P{t)ipt at later times. The adiabatic theorem as- 
serts however that such a property holds approximately 
in the sense that 



- P{t))M = 0{e) 



(2) 



with ll'^ll = {tpl^y/^ the length of the vector ip. 

While the error estimate (||) is undoubtedly correct, it 
really begs the question, since the nature of 0{e) is left 
unspecified. It could be that a piece of size e of V't leaks 



out into the orthogonal subspace (1 — P{t))H. Alterna- 
tively \ipt){flJt\ is slightly tilted relative to P{t). As first 
recognized by Lenard ^ , and on a more refined level in 
lUi H' 0' §i it is the latter option which is realized by 
the solution to (|l|). There is then an iterative procedure 
for constructing a projection £-close to P{t), such 

that ipt remains in P^{t)TL up to any order, in symbols 

\\{l-Pm))M=0{en, ^'o^P'mo. (3) 

The power series for P^{t) is asymptotic, with an error 
of order exp[— 1/e] in case H{t) is analytic in t p^ . 

Such results are of physical interest only if higher or- 
der corrections can be computed in a systematic fash- 
ion. We explain a scheme which naturally generalizes 
to the space-adiabatic decoupling. The idea is to sep- 
arate kinematics, the subspace P^ [1)1-1, from dynamics, 
the unitary evolution inside P'^{t)Ti as generated by some 
effective Hamiltonian. For this purpose we choose a ref- 
erence subspace of dimension m with time- independent 
basis Ixq,) , a = 1, m. The effective Hamiltonian oper- 
ates in this subspace and to order is given by 



SapE{t)-ie{ipa{t)\ipp{t)) (4) 
{l-P{t))Mt)) + Oie^), 



+ \e^^^mH{t)-E{t))-^ 

a,(3 = l,...,m. Note that the dynamical problem (|^) 
has been reduced to the dimension of the relevant sub- 
space. Secondly there is a unitary, U^{t), which rotates 
the reference subspace into P^{t)'H. To order e one finds 

a=l 

+ze{H{t) - Eit))-\1 - Pm^^it))) (x„| + Oie') . 

Thus, if x'^i't) denotes the time-evolved vector in the 
reference subspace, i.e. the solution to iedx'^{t)/dt — 
H''^{t)x^{t) , x''(0) = X, then the 0(£°°)-approximate so- 
lution to (m) reads 



(6) 
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which is vaHd on the slow time scale with an error 0{e'^) 
upon inserting (^) and 

Physically, the separation into prescribed slow degrees 
of freedom and rapidly adjusting fast degrees of freedom 
is somewhat artificial. In actual fact there is a single 
time-independent Hamiltonian governing all degrees of 
freedom. For example, the molecular Hamiltonian is in- 
dependent of time and the separation into slow nuclei and 
fast electrons results from the small mass ratio. What is 
then needed is a spatial version of the time-adiabatic the- 
orem in the sense that some dynamical degrees of freedom 
are frozen out. While the physical mechanism behind 
such a space-adiabatic decoupling may differ widely, our 
discovery is that there is a unifying and general theoret- 
ical framework. We like to compare the situation with 
standard bound state perturbation theory. It is always 
the same algebra, but the physical manifestations are rich 
and differ widely depending on the context. The general 
scheme to be sketched below is applicable to a variety 
of systems such as the Dirac equation with slowly vary- 
ing external fields, the Born-Oppenheimer Hamiltonian, 
electron motion in periodic potentials, magnetic Bloch 
bands, motion on embedded manifolds with constrain- 
ing transverse potential, and others. These systems have 
been studied before by model dependent methods and 
it is not uncommon to have a debate as regards to the 
structure of second order corrections |ll, 12, p^. 
In contrast our perturbation scheme systematically com- 
putes higher order corrections and provides error bounds. 
Such a scheme was envisioned by E. I. Blount in a remark- 
able series of papers Particular aspects are 
studied by Littlejohn and Flynn ^ and Nenciu and 
Sordoni The full justification of the space-adiabatic 
perturbation theory is beyond the size of a letter and 
given elsewhere [^^. Here we merely explain the struc- 
ture of the leading order terms. As novel applications we 
discuss the second order Born-Oppenheimer theory and 
derive the gyromagnetic ratio (g-factor) for an electron 
coupled to the quantized radiation field within nonrela- 
tivistic quantum electrodynamics. 

We claim that the spatial version of Eq. Im) reads 



(7) 



q G IR'' is a position like variable and p the corresponding 
canonical momentum. Again, e <C 1 is a dimensionless 
scale parameter, which controls the variation of q and 
p. The phase space for the slow degrees of freedom is 
r — H{q,p) is an operator valued function on P, 
which acts on the Hilbert space H{ of "internal" degrees 
of freedom, and q = x,p = — iV^. To properly define 
H{q,ep) we adopt the Wcyl rule for operator ordering, 

H{q,ep)i;{x) = {2tt)-'^ J dVd'*Ce*(^-^')-« 

xH{^{x + x'),e0^ix'), (8) 



where H acts as an operator on ip. The "spinor" wave 
function ipi^) takes values in Tif. Thus the Hilbert space 
for the Schrodinger equation (0) is H = L'^(R'^,n{) = 

The eigenvalues of the Hamiltonian H{q,p) define the 
energy bands. One of them, or one group of them, is 
considered to be physically relevant and, only for the 
sake of simplified presentation, is assumed to consist 
of a single band energy E{q,p) of constant multiplicity 
m. Thus there are m orthogonal eigenvectors satisfy- 
ing H{q,p)\ipa{q,p)) = E{q,p)\ipa{q,p))- They span the 
projector P{q,p) = jy^=i \'>Pa{q,p)){'^a{q,p)\- As crucial 
assumption one needs that for each q,p the energy E(q,p) 
is separated by a gap from the remainder of the spectrum 
of H(q,p). In spirit, to lowest order in e, the projector 
onto the adiabatically decoupled subspace of Ti. is the 
Weyl quantization P{q,ep). However, P{q,£p) is not a 
projector, in general. Still on abstract grounds one can 
construct a time-independent subspace P'^'H, the projec- 
tion being e-close to P{q,ep), which is adiabatically 
decoupled to all orders in e, in the sense that 

- P')M ^ 0{e^) , yJo = P'i'o, (9) 

for tpt the solution to (0), compare with (^. 

To systematically compute approximate solutions we 
disentangle kinematics and dynamics. The reference 
Hilbert space is now L^{R )(g,C"' with fixed spinor basis 
\Xa),C( = 1, ...,m. The effective Hamiltonian to order e 
turns out to be the Weyl quantization of 

Hfpiq,p) = S^f,Eiq,p) (10) 
-^ls{MQ,P)\{Hiq,p) + E{q,p),i^p{q,p)}) + 0{e^) . 

{•, •} is the Poisson bracket defined for the operator- resp. 
vector-valued functions A, ^ as {A, ip} = VpA ■ Vqip — 
\7qA-Vpip. The reference subspace is rotated by the uni- 
tary If^ into the adiabatically decoupled subspace P^Ti. 
To lowest order one obtains the Weyl quantization of 



Uo{q,p) = ^ \i'a{q,p)){xa\ 



(11) 



Uo{q,ep) is not unitary in general, but the error is in 
higher orders. As in the time-adiabatic setting the ap- 
proximate solution to the Schrodinger equation (|^) is 
given by 

= C/^(g,£p)exp[-zi/^*f(g,£p)i/£]xo + , (12) 



where ipQ = U^{q,ep)xo- Inserting from ( |10[ ) and (|l] 
the solution holds on the slow time scale with an error 

Higher order corrections to (|o|) , (^l]) are available , 
but they rapidly increase in complexity. To provide some 
idea we explain the second order Born-Oppenheimer ap- 
proximation, for which the Hamiltonian has the form 



HBo{q,p) 



V{q), 



(13) 
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where 1 is the unit operator and V{q) is some opera- 
tor both acting on Tif. For a molecule, q would stand 
for the coordinates of the nuclei and V{q) is the elec- 
tronic Hamiltonian with the nuclei fixed at q. The adia- 
baticity parameter is e = (wei/mnuc)"'^^^, assuming that 
the nuclei have all the same mass, and Hbo(J1j£p) = 
— -I- V{x) is the standard molecular Hamiltonian. 
The relevant subspace is spanned by eigenvectors to 
Viq), ViqMUQ)) = e(g)|V'a(g)),a - l,-,"!. Thus 
E{q,p) = ^p^ + e{q) and P{q,p) depends only on q 
as P{q) = \'^a{q)){ipa{q)\, which implies that also 

P{q) is a projector without additional corrections from 
the operator ordering. Generalizing jl^, ^ , the effective 
Hamiltonian including second order is computed to 

ni 
7 = 1 

+ e{q)Sa,p + ^e^{Vi>o.{q)\-V^Pp{q)) 

- e^{p • V^„(<z)|(V^(g) ™ eiq))-\l - P{q))p ■ V^PM) 

+ 0(£3), (14) 

where the "vector potential" of the Berry connection 
reads Aap{q) = i{^a{q)\'^'ipp{q)) ■ As before, the effective 
quantum Hamiltonian to that order is the Weyl quanti- 
zation H^{q,ep). For the unitary, which rotates the 
reference space back to P'^TC, one obtains 

m 

UBo{q,p) = {\Mq)) + i£{V{q) - e{q)y' 

a=l 

X (1 - P{q))\p ■ ^Mq))) iXcl + 0{e') . (15) 

Note the similarity with the timc-adiabatic approxima- 
tion d), d). 

We skipped over the issue on which time scale ( p^ is 
valid. Clearly, there are two different ways of reading 
our result. Either one considers the time scale 0(e°) 
on which (^ is written. Then, inserting ( p^ and (15), 
the error in ( p^ ) is 0{e'^). Or one considers the longer 
time scale ©(e"^). Then precision is lost and the error 
is guaranteed to be 0{e) only. Note that in this second 
variant the adiabatic decoupling time is much longer than 
the one which would result from a semiclassical analysis 
of the Hamiltonian (g, ep) . Its validity is restricted 
to the Ehrenfest time which usually is of order log£~^ in 
our units. 

To have a concrete physical example we discuss the gy- 
romagnetic ratio of the electron in nonrelativistic QED. 
The electron is subject to a weak uniform external mag- 
netic field. The internal degrees of freedom are the pho- 
tons which rapidly adjust to a state of lowest energy con- 
sistent with the momentary state of the electron. The 
orbital motion is circular with frequency ujc and the spin 
motion is precession with frequency LOg. As the standard 
definition of the g-factor we adopt 



To apply our general scheme it is convenient to con- 
sider first the motion in arbitrary external potentials of 
slow variation. The (7-factor will be retrieved as a spe- 
cial case. In dimensionless Heaviside-Lorentz units the 
Hamiltonian with zero external potentials reads 

H = il/2){a ■ (pci - eA^{x)))^ + Hf , (17) 

a = e?/Aii the fine-structure constant |Q. a;,Pci are 
position and momentum of the electron, a is the vector 
of the Pauli matrices, Hi — Ea=i / d?'k\k\a^ [k^ X)a{k, A) 
is the energy of the photons with dispersion uj{k) — \k\ 
and helicity A, and A^p[x) is the quantized transverse 
vector potential, which is smeared with the form factor tp 
(=ultraviolet cutoff). The dressed electron responds to 
the slowly varying external potentials (pcxisx) ^ Acy,{ex) 
through a slow variation of the total momentum p — 
Pol + -Ff 7 where Pf = Ea=i 2 / (Pkka^{k, X)a{k, A) is the 
momentum of the photons. To write the Hamiltonian in 
the standard space-adiabatic form (^ one thus has to 
substitute p — Pf for poi and x/s for x in order to arrive 
at the operator valued function 



H{q,p) = i(p-Pf-eA^-e 



.^ex(g)) (18) 

^e{Bp + eBcx(<?)) • a- + e0ox(g) + Hf 



A 



A^{0), B, 



'■p 



X 



with the abbreviations 
AM)),Bc^{q) = V X Acy^{q). The Weyl quantization of 
(O) is the Hamiltonian of nonrelativistic QED. H{q,p) 
acts on the electron spinor and the wave function for the 
photons. Thus Hf = ® J- with T the photon Fock 
space. Note that H{q,p) has a rather particular struc- 
ture which we abbreviate as 

H{q, p)=D{p)- ieePex(g) • <J + ec^^q) , (19) 

where p = p — eAc-^{q) is the kinetic momentum. 

The physically relevant subspace is spanned by 
the eigenvectors of lowest energy. They satisfy 
■D(p)|V'q(p)) = ^^o(p)IV'a(p))> which implies E{q,p) = 
Eo{p - eylox(g)) + e(/)ox(g)- By spin degeneracy the 
ground state subspace turns out to have dimension 2, 

In addition, Eq(j)) is isolated 



I.e. m 



,a 



±1 



25| 



2Ws/Wc 



(16) 



from the rest of the spectrum for \p\ smaller than Comp- 
ton momentum, provided the interaction is cut off in the 
infrared. We want to fix a basis such that V'+ corresponds 
to the electron spin pointing along the +z direction. For 
this purpose one notes that the total angular momentum 
J in the direction of p is conserved, [D{p),p ■ J] = 0. 
Thus the eigenvectors along ±p are uniquely fixed and 
by suitable linear combination therefore also along ±z. 

At this stage one merely inserts in (|lo|). The effective 
Hamiltonian including first order is computed to 

Hfpil'P) = (^^o(p-e^ex(g)) + e0cx(g))(5a/3 
-^eeBcxiq) ■ (V^a(p)|cr'0/3(p)) 
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-ee( - V(/)cx('?) + v X Sox(g)) ■ (V'q(p)I«VV'/3(pO) 



(20) 



where (-I-) is the scalar product in C^®^, p = p— e^ex(<z) 
as before, u = \/Eq{p). Upon Weyl quantization iJg^ 
is the Peierls substitution. Hf^ describes the spin pre- 
cession. The effective Hamiltonian ( pO| ) is most conve- 
niently studied through semiclassical methods |2^]. The 
center qt , pt of the wave packet is governed by the clas- 
sical Hamiltonian Eq{p — eAcx{q)) + e(f>c->L{q)- With our 
choice of time scale, qt,Pt no longer depend on e. In- 
serting a particular orbit qt , pt in H^^ and cancelling the 
e-factors in (0) and @) yields the spin precession as 



*^Xt = Hf{qt,pt)xt ■ 



(21) 



Orbital motion and spin precession vary on the same time 
scale and the spin does not react back on the orbit to our 
order. 

Of particular interest is a uniform external magnetic 
field, Bcy^iq) = -B, (?I)cx('z) = 0, and small velocities v = 
\7Eo{p), which means pt = in good approximation. 
Let us introduce the bare mass m of the electron, the 
velocity of light c, and Planck's constant h. Then the 
spin Hamiltonian becomes 



rrspi 

^af3 



-^B.[{MO)\<^MO)) (22) 



which implicitly defines the frequency ojs — ^AB\ of the 
spin precession. On the other hand the orbital motion is 



rricsq ^ -q x B , 
c 



m,tf = 1/K(0) , 



(23) 



which yields the cyclotron frequency ujc — e\B\/mcf[C. 
Together with the definition (|l6|) one deduces the non- 
perturbative expression of the g-factor valid for any value 
of the fine structure constant. 

The matrix elements in ( |2^ ) and for Eq{0) are not 
available as closed formulae. The only possibility is to 
expand in e around e = 0. To second order in e and 
removing the ultraviolet cutoff one obtains 



(24) 



in agreement with pq |, who use static second order per- 
turbation for the Zeeman splitting. 

In summary, the expressions (|l^) and (|l^) and their 
higher order corrections provide a systematic approxima- 
tion in space-adiabatic situations. To compute them one 



needs as explicit as possible the eigenbasis in the rele- 
vant subspace. Just like in static perturbation theory for 
eigenvalues at second order, also here the inverse oper- 
ator {H{q,p) — E{q,p))~^ appears, compare with (|l4|), 
(p^), which is not necessarily readily available. From our 
own experience approximations beyond second order be- 
come tedious and most likely carry no additional physical 
insights. 
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